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ABSTRACT 


; The LLNL Paladin FEL experiment is shown to exhibit clear and dramatic 
effects governed by the electron beam velocity distribution for the first time. The 
FEL integral equation is used to show that there is significant broadening of the 
gain spectrum due to the Gaussian velocity distribution, and also shows a plateau 
in the gain evolution along the undulator due to a triangular-shaped velocity 
distribution. 

The gain spectra and power evolution from simple, single-mode simulations 
are compared to the LLNL ELF experiments. The microwave power evolution 
along the undulator is compared as well for both the tapered and untapered 


undulators. In all cases, the agreement is found to be good. 
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I. INTRODUCTION 


In a Free Electron Laser (FEL), a relativistic high-quality beam of 
electrons, passing through an undulator magnet, generates a powerful and 
coherent electromagnetic wave [1]. Coherence may develop after many 
passages in the undulator, when the device operates as an oscillator, or in 
one passage in the amplifier configuration. The choice between an oscillator 
configuration and a amplifier configuration is generally determined by the 
accelerator that drives the FEL; specifically by the pulse format and peak 
current of the accelerator. Induction linear accelerators produce short 
pulses, usually 10-100 ns long, that are widely spread in time with high 
peak current. The pulses are too far apart to be synchronized with optical 
pulses in a resonator. To be interesting, the amplifier must have high gain 
in a single pass of the light through the amplifier. The high peak current of 
an induction linac permits both high gain and high efficiency, the fraction of 


electron beam energy converted to light. 


The evolution of coherent electron bunching is the key element to any 
FEL gain mechanism. In a high-gain amplifier, like the Paladin FEL [2-4] 
at Lawrence Livermore National Laboratory (LLNL), maintaining coherence 
of the electron bunch over a significant interaction length imposes important 
restrictions on the initial electron beam quality. The design of an high gain 
FEL amplifier requires high beam current and a long undulator; both 


features tend to allow gain to be easily degraded by poor beam quality. 





Paladin is the first short-wavelength FEL to operate in the high current 
regime where the growth mechanism is dominated by the longitudinal 
velocity distribution of the electron beam. It is observed that the growth of 
the optical power along the 15 m Paladin undulator is not exactly 
exponential, but is determined by the shape of electron distribution 
function. Similarly, the gain spectrum measured in the earlier 5 m long 
Paladin amplifier shows a dramatically broader and lower shape due to the 
electron beam distribution. The FEL integral equation [5-7] is used to 
calculate the optical field evolution in the presence of an arbitrary electron 


velocity distribution. 


The LLNL Electron Laser Facility (ELF) FEL is a single-pass, high- 
gain, high-power amplifier operating in the microwave regime at 34.6 GHz 
frequency [8-11]. The characteristics of this FEL have been extensively 
studied with large simulations including many transverse modes, 
diffraction, beam emittance, and betatron focusing. A relatively simple, 
single-mode model of the FEL interaction using only the electron pendulum 
and optical wave equations can describe the most important features of the 
ELF FEL. The waveguide, used in the experiment to confine the natural 
diffractive spreading of the microwave beam, can be simply modeled by a 
filling factor and a shift in resonance. 

The most significant contributions to the analysis of a high-gain, high- 
power FEL are summarized below. 

° Gain spectrum broadening of the 5 m Paladin experiment due to 


poor beam quality are examined by FEL integral equation. 
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Plateau in the 15 m Paladin FEL gain evolution is explored by 
the velocity distribution functions. 

The asymmetrical shape of high-gain ELF FEL is examined by 
simple FEL equations. 

Gain evolutions with both untapered and tapered ELF FEL are 


explored using simple FEL equations. 








II. BASIC FEL THEORY 


A. ELECTRONS IN THE UNDULATOR 


In order to discuss the principles of the FEL, it is important to consider 
the individual electron trajectories in the undulator field alone, without the 
presence of light. The interaction between the negatively charged electrons 
is negligible in typical FELs [1], and the beam is focused externally before 
entering the undulator. Under the influence of the undulator magnetic 
field, the electrons will deflect periodically in the transverse direction, or 
"wiggle" as they travel along the optical mode axis. The ideal form of the 


on-axis helical undulator field atx =y =0 is 
B = Blcos(kz),sin(koz),0]_, for O<z<L , (1) 


where B is the peak undulator field strength, kg = 2n/A, is the undulator 
wavenumber associated with the wavelength Ao, z is the distance along the 
longitudinal axis of the undulator, L = NA» is the undulator length, and N 
is the number of undulator periods. The Lorentz force equations [12] 


become 


i oo oe ae 
at) = ——(BxB) , and —F=0, (2) 


where y= (1-B 2)-12 is the Lorentz factor, m is the electron rest mass, 


e =lel is the electron charge magnitude, c is the speed of light, 








B=v%c = BB, B, is the perpendicular component which refers to the 
transverse direction, and B is the parallel component along the axis of 
undulator. We see that the electron energy, yc”, is a constant of motion in 


the magnetic field above. Substituting (1) into (2) gives 
d meee Ae hed 
By =e B, B[-sin (k oz ),cos (Rk oz ),0) (3) 


The velocity can be solved exactly from (3), but assuming perfect injection, 


y > 1, and B, = 1, the velocity components can be written as 


Ba = Een (k oct )-Asin (k oct Be : (4a) 
1K? 1/2 
where 8, = pe) ; (4b) 


and K =eBi /2nmc? is the dimensionless undulator parameter or vector 
potential. 
The electron position can be found by an additional integration of (4a) 


and is 





Kho 

Oix 

t),—— t),B,ct| , 
oe ye oy COS (Wot ), Bc 


Kh 
Pt) = | 


where @) = B,koc. In a typical FEL, K = 1, A, = 5 cm, andy = 100, so that 
the electrons travel at speed B,c =c along the z-axis for several meters 
while executing small transverse oscillations with amplitude KA)/2zy = 0.1 
mm, which is smaller than the typical beam radius. The transverse 
deflection of the highly relativistic electron beam in the alternating, periodic 
field of the undulator leads to spontaneous emission. Even if the deflection 
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is too small to be visualized on the scale of the undulator, it is still 


responsible for crucial spontaneous emission, and FEL gain. 


The spontaneous emission mechanism is essentially classical, and can 
be described by the Larmor formula giving the total power emitted from a 
single electron as P, = 2e >A 7/3c = 2yB*cr2/3, where r, =e7/mc? is the 
classical electron radius (in cgs units) [12]. This radiation is emitted into a 
small cone about the forward direction with angular width = y!. This 
narrow-spectrum radiation is trapped in the optical resonator and used to 
supply the starting photons for oscillator operation. The interaction 
between these photons with the electron beam in the undulator is the 
essential concept of FEL physics. In the amplifier configuration, the 


photons are supplied by the seed laser. 


The fundamental feature of the FEL mechanism is that the beam 
electrons undergo axial bunching in the combined undulator and radiation 
fields to exchange energy with an optical field. Assuming a circularly 


polarized plane wave, the optical fields are taken to be [12] 
E. = E(cosy,-siny,0) , and B. =EH(siny,cosy,0) , (5) 


where E is the optical field strength, y = kz-wt+ with carrier frequency 
w=kc, k = 2n/d is the optical wavenumber correlated with the wavelength 
A, and 9 is the optical phase. 

Since the electron travels at slightly less than the speed of light, several 
optical wavelengths pass over an electron as it travels the full length of the 


undulator. When exactly one optical wavelength passes over an electron as 


that electron travels through one unduiator wavelength, resonance is 
established. Efficient energy exchange requires that the electrons 
experience near resorant forces from the optical and undulator fields. For 
y 1, this resonance is achieved when the optical wavelength, the beam 
energy, the undulator strength, and the undulator period approximately 


satisfy [1] 


(14K) 


27 


The expression for the resonant FEL wavelength above shows one of the 


A= 


FEL’s most important attributes, continuous tunability. As the electron 
energy from the accelerator is varied, or as the undulator strength K is 


varied, the resonant wavelength is changed. 


B. OPTICAL FIELD EVOLUTION AND THE WAVE EQUATION 
In order to follow the growth of the optical field it is necessary to 
develop a wave equation. Start with Maxwell’s wave equation in the 


Coulomb gauge [12], 





where A’ is the optical vector potential, and J , is the transverse current 
density. An appropriate form of the vector potential for a circularly 
polarized optical wave is A =(E/k [sin y,cos y,0]. The total current can be 
represented as the sum of the single-particle currents in the electron beam, 


i = -€C YB SPP ), where ® follows the trajectory of the i‘ particle. 
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Assuming that the amplitude and phase of the optical field is taken to vary’ 
slowly over an optical period (E <a, 6 <@o) and over an optical 
wavelength (E’ «<kE,o’«<ko) and keeping only terms with single 


derivatives, the equation can be written as 


[2 + 12 |pte = ta etS> : where C= (Rotk )z — ot 
oz cat y 


p is the electron beam particle density, the brackets, <..>, indicate a phase 
average over the sample electrons within a wavelength of light, and € is a 
dimensionless variable that describes the electron’s phase with respect to 
the combined optical and undulator fields within roughly one optical 
wavelength, since k »k, and z =8,ct. If we assume a long electron pulse 
and don’t consider longitudinal optical modes, we can remove z -dependence 


from the wave equation, then the wave equation can be written as 

° . 

a=-j<e'S> , for O<t<l1 , (6) 
where the complex dimensionless field envelope and current density are 


g = 40eKLNEe!® j = SN(exKL)’pF 
mc? ymc2 
F =r,2/w? is the "filling factor" which is the ratio of the electron beam area, 


2 


nr2, to the optical mode area nw’, averaged over the undulator length, 


° 
t=ct/L is the dimensionless time, and the " " denotes a time derivative 
with respect to t. If many electrons are randomly spread in longitudinal 


position, then the phase average < e?$>=0, with no optical field growth. 


In order to understand the FEL interaction, first we must examine how 


electrons evolve from their initial positions within a wavelength of light. 


C. ELECTRON DYNAMICS AND THE PENDULUM EQUATION 


The interaction of the electrons and the optical fields inside the 
undulator is described by the complete Lorentz force law and the energy 


equation [12]; 


dm. e 
08) = B+ BB) , (7a) 
and: 2h ee ee (7b) 
t mc 


Inserting the full combined fields of the undulator field (1) and optical field 


(5), the transverse motion with perfect injection can be found to be 
B, =- 7 [eos (R oz ),sin (& oz ),0) + “isin y,cosy,0) , (8) 


where A =eER2nmc?2 is related to the optical vector potential and is 


analogous to the definition of K. 
From (7b) and (8), the change of electron’s energy is found to be [1] 


pee costae). (9) 
a 


c 





The exchange of energy of the electrons is easily understood in terms of the 
electron’s phase, €. For example, electrons with phases such that, 
—n/2 < (C+) < n/2, will get energy from the optical field causing them to 
accelerate. Other electrons with phases, 1/2 < (C+) < 3x/2, will lose energy 
to the optical wave causing it to grow. For a random distribution of initial 
phases, about half the electrons gain energy and about half lose. But, under 
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the influence of the combined fields, the electrons will begin to "bunch" at 


desirable phases and cause gain during a pass down the undulator. 


To further clarify the dynamics of the electrons in the undulator, it is 
convenient to rewrite (9) in terms of € alone. Using the relation 
y?= 1-B?-B2, and the conditions y» 1, B, = 1, and B? = Kr, it can be 
shown that the individual electrons follow motion described by a simple 


pendulum equation [1] 


tole) ° 


C =ve=lalcos(C+o) , (10) 
where v = t = L[(k +k 9)B, -k ] is the dimensionless phase velocity. 


The wave equation (6) and the pendulum equation (10) are valid in both 
strong and weak optical fields with high or low gain. They are not accurate 
for the high-efficiency FEL where there are significant changes in the 
electron energy ync”. The electron energy is treated as a constant along the 
undulator in the definition of the dimensionless field strength la |, and is not 
followed self-consistently. The optical field strength la! explains why 
stronger optical fields produce rapid electron evolution and, how more 
relativistic electrons, with a larger y, require a stronger optical field E to 
accomplish the same bunching. When an electron loses energy, the value of 
la | increases so that the coupling between ¢ and v increases. When more 
than 10% of the electron beam energy is converted to radiation, the FEL 
dynamics can be altered so that the more complete wave and pendulum 


equations may be used. 
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III. BEAM QUALITY AND GAIN EVOLUTION 


A. THE FEL INTEGRAL EQUATION 

To describe the FEL classically, the electron dynamics are governed by 
the optical wave (6) and electron pendulum (10) equations. In the case of 
weak fields, lal «<a, FEL dynamics can be expressed in the form of an 
integral equation by solving the electron pendulum and optical wave 


equations self-consistently [1, 5-7], 
7 iy . o gl , -I Vot’ , 
a(t) = on Jat vF(re a(t-t’) , (11) 


where F(t’) = fdaf (q)e#* is the characteristic function of the distribution 
f(q), f(q) is the distribution of initial electron phase velocities v; = vo+q 
about vo, and [f(q)dq =1. Note that F(t’) and f(q) only depend on the 
electron beam’s initial conditions with subsequent dynamics described 
exactly. The FEL integral equation (11) is generally solved numerically, and 
has a distinct advantage over simulations which must follow many tens of 
thousands of sampled electrons. 

When the current is low, j < x, the field a (t-t’) can be approximated by 
the initial field a(0) =a g in the integrand of (11). In the case of perfect 
beam quality, f (¢) = 6(0) and F(t’) = 1, the field at the end of the undulator 


a(1) is easily evaluated to find the spectrum for low-gain, 
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G(vo) = j [2 — 2cos (vo) — vosin(vo)] / ve (12) 


where gain is defined as G = (a%(1) - a2) /a@. The antisymmetric shape of 


the low-gain spectrum has width Avy = 2x and a peak value G = 0.13) at 
Vo = 2.6. These features are among the first FEL characteristics found 
experimentally [1]. 

When the current is high, j » 2, the optical field grows exponentially 
and at resonance, Vo = 0, where the growth rate is maximum, the gain is 


G(t) = exp ((j/2)7V3t) / 9. Near resonance, Vo = 0, the final gain at t = 1is 
G (vo) = exp (V3(y/2)"3) (1 — vg/3V8Uj/2)"3) 79. (13) 


The gain spectrum for high current is symmetric about resonance, and 
drops off at a_ characteristic value of the phase _ velocity 
Avo = 2(3V3)"2(j/2)"6 = 4.2276 For large j, say j = 30, the peak gain, 
G = exp(V3(j/2)"/3y/9 = 8 is much larger over a broader spectrum, AVo = 7.5, 


than in the low-gain case. 


B. ELECTRON DISTRIBUTION FUNCTIONS 


Many experiments do not have the perfect beam quality used in the 
analytic results (12) and (13). Depending on specific experimental 
conditions, the actual shape of the electron beam distribution function may 
vary from one experiment to another. An initial electron phase velocity 
v; =L{(k +k o)B,(0)-k] depends on the initial electron z velocity cB, (0), 
and is therefore a function of both the energy spread and angular spread in 


the beam. Consider a relativistic electron beam with average energy ymc2, 
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and spread described by width og = 4nNAy/y. The resulting distribution to 
be used in the integral equation (11) is the Gaussian distribution [5-7], 


e -q7/20% 


oa. 
f(q) Vinag or all q (14) 


Alternatively, an angular spread caused by emittance in the beam can 
result in the exponential distribution for f(q). <A relativistic electron 
entering the undulator with velocity cB) at a small angle 6 has a reduced z 
velocity, cB, = cBycos® = cBo(1-92/2), and the corresponding change in the 
phase velocity, v(0) = vo — 2nN y°02/(1+K?). A Gaussian spread in angles 0 of 


width A@ results in the exponential distribution [5-7], 


~q /9 


_ @ 
f(q)= Bs 





for q <0 and f(q)=0 for q>0O , (15) 
) 


where 6g = 4nN y7A@2(1+K“). If a very broad angular distribution is cutoff 


by an emittance filter, the resulting phase velocity spread is a triangular- 


shaped distribution with positive slope [7], 


f(q)= (0, +9) for -o, <q <0 and f(q)=0 otherwise , (16) 
0 


A 
where o, is the width of the triangle. 

In a realistic FEL, the correct shape of f(q) is generally not known 
precisely, and is determined by several complicating factors involving the 
accelerator and electron beam transport properties. The L = 5m Paladin 
experiments are taken to be described by the Gaussian distribution (14), 


whereas the L = 15m Paladin experiments made use of an emittance filter 
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and transport system that resulted in a distribution more like (16) or 


possibly (15). 


C. GAIN SPECTRUM BROADENING IN THE 5 m PALADIN 


As the LLNL Paladin experiment developed, measurements were made 
with undulators of first 5 m and then later 15 m length [2-4]. In the 5 m 
experiment, the electron beam current was about 600 A at about ymc? = 41 
MeV energy and had a radius of r, = 5 mm. The L = Nay = 5 m undulator 
contained N = 62 periods of wavelength A, = 8 cm, and undulator parameter 
K =0.9. The radiation wavelength was 1 = 10.6 um in a radiation mode 
with waist wo = 4 mm that was slightly smaller than the electron beam 
radius r, = 5 mm at the beginning of the undulator so that it reduced the 
effective current density. Natural diffraction allows the optical mode to 
expand, and we determine the filling factor F in J by averaging along the 
undulator. As a result, we describe the 5m Paladin experiment by 
dimensionless current j = 30. The observed electron beam brightness 
establishes that the spread in initial phase velocities is about Avg = 12, but 
the shape of the distribution is not known. We take the shape to be the 


Gaussian distribution as described in (14) with og = 12. 


In figure 1, the weak-field gain spectrum (shown as dots) was measured 
by varying the undulator field strength B over a small range to change the 
resonance condition Vy [5-7]. The observed gain spectrum width is Avg = 40, 
with peak gain of G = 12%. The solid curve uses the integral equation to 


evaluate gain with j = 30 and 6g = 12. The theoretical gain spectrum was 
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observed to be antisymmetric in shape with a peak gain of G = 0.15. The 
observed and theoretical gain spectra are dramatically distorted by the 
effects of poor beam quality. Figure 2 shows, for comparison, the gain 
spectrum determined by the integral equation with j = 30 and og = 0 
describing perfect beam quality. Poor beam quality is responsible for 
decreasing the FEL peak gain from G = 9 in figure 2 to G = 15% in figure 1, 
and distorting the gain spectrum from an approximately symmetric shape of 
width Avy = 6 to an antisymmetric shape of width Avg = 70. Note that the 
observed and theoretical antisymmetric spectra in figure 1 are both much 
broader than the antisymmetric low-gain spectral shape in (12) with width 
AVo = 2n. 

The shape, overall width, and peak value of the 5m Paladin gain 
spectrum is dramatically altered from the idealized gain spectrum (either 
(12) or in figure 2), and in this case, is determined by the initial distribution 


of electron velocities (14). 
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Figure 1. Gain evolution of the 5 m Paladin with poor beam quality 














Figure 2. Gain evolution of the 5 m Paladin with ideal beam quality 
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D. A PLATEAU IN THE GAIN EVOLUTION OF THE 15 m PALADIN 


We also identify characteristics of the Paladin 15m experiments that 
are unique to the initial electron distribution function shape. An emittance 
filter was used to improve the beam quality entering the undulator, and a 
detailed study of the filter and transport system concluded [13] that the 
distribution function shape should look something like the exponential or 
triangular distributions, (15) or (16). 

In the 15 m experiment, the electron beam current was about 500 A at 
ymc2= 45 MeV energy. The longer undulator was L =15 m _ long, 
containing N = 187 periods with K = 1.0. The high gain observed in the 15 
m experiment causes the optical mode and electron beam to be well matched 
over the whole undulator length, so that we take F = 1. Other parameters 
are the same as in the 5 m experiment. Because of the longer undulator 
length, the 15 m experiment is described by large j = 1400, and the 
triangular distribution (16) with o, = 13. For such a large current, the ideal 
gain, with perfect beam quality, o, = 0, would be G = 5x10° peaked at vg = 0 
and would evolve exponentially in t. Figure 3 shows the theoretical gain 
evolution, G(t), found by solving the integral equation with the triangular- 
shaped distribution (16) of width o, = 13, 7 = 1400, and vy = 0. The gain 
evolution, G(t), represents the optical power growth along the undulator 
length as t goes from 0 to 1, and it shows a dramatic plateau about half way 
along the undulator from t= 0.5 to 0.7. The choice of an exponential 
distribution (15) shows significantly smaller gain than the triangular-shaped 


distribution and with a more subtle plateau. 
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Figure 4 shows the optical power, P, measured along the 15 m 
undulator length [2], corresponding to t=0—71. The changes in power 
indicate the gain evolution, and the figure shows a the prominent plateau in 
the middle of interaction length. It appears that the shape of the 
distribution function in the Paladin 15 m experiment has resulted in a 
plateau in the FEL interaction half way through the undulator, and the 
experimental plateau in the gain evolution is verified by the integral 
equation with the positive sloped triangular-shaped distribution. 

In conclusion, the results of the Paladin 5 m and 15 m experiments, 
have provided a unique opportunity to observe how the FEL gain spectrum 
and gain evolution can be determined by the electron beam velocity 
distribution instead of the natural high-gain mechanism. The 
measurements provide the first comparison between theory and experiment 
in a short-wavelength device where the electron beam distribution plays a 
dominant role in the FEL mechanism. The gain spectrum shape, width, 
and evolution can be used to determine unique information about beam 


quality. 
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Figure 3. Theoretical gain along the 15 m Paladin undulator 


20 















Power (W) 


Figure 4. Experimental gain along the 15 m Paladin undulator 





IV. ASIMPLE MODEL OF THE ELF I EL 


A. DESCRIPTION OF ELF EXPERIMENT 


The LLNL Electron Laser Facility (ELF) FEL is an amplifier with very 
high gain operating in the microwave regime using the Experimental Test 
Accelerator (ETA). The experiment examined high-power saturation, and 
the effect of tapering the undulator field. ETA uses a 6 kA electron beam 
with a normalized emittance of 1.5n cm-rad [8]. An emittance filter reduces 
the beam current to approximately 450 A over a pulse length of 15 ns with 
normalized edge emittance of 0.47n cm-rad. The beam energy is 3.3 MeV 


3 


with y = 7.5. The electron beam particle density is p = 2 x 10!! cm™, and 


the beam radius is r, = 4 mm. 


The L = 3m long ELF undulator is linearly polarized, with wavelength 
Ay = 9.8 cm and undulator constant K = 2.8. The peak field strength is 
B = 4300 G over N = 30 undulator periods. The Rayleigh length is less 
than 1% of the undulator length so that ELF must use a 3 cm x 10 cm 
rectangular waveguide to confine the radiation wavefront around the 
electron beam for continued amplification. Without a waveguide, natural 
diffraction would spread the optical wave cross-section away from the 
amplifying electron beam and reduce coupling. The waveguide cutoff 


frequency is 5 GHz in the TE», mode, so that the FEL at 35 GHz operates 
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far above cutoff. 





B. SIMPLE FEL EQUATIONS FOR ELF EXPERIMENT 


The FEL equations described in previous chapters are the self- 
consistent pendulum and wave equations. They provide a simple model of 
the FEL that is valid in strong and weak optical fields for high and low 
beam current. The assumptions required to derive the simple form of the 
equations are y> 1 and N »1. Since the most important aspects of the 
interaction take place near resonance, the changes in the electron energy, 
Ay/y, are small when N >1, and the FEL efficiency is low. For high 
efficiency, in the tapered undulator or in a short undulator with few periods, 
the changes in the electron energy are expected to be large. High efficiency 
FELs the undulator can be tapered to maintain resonance with the trapped 
electrons that lose energy to the optical wave. The untrapped electrons will 
drift far off-resonance. Tapering and high-efficiency can be added to the 


basic FEL electron equation of motion (10), and it becomes [14, 15] 
© =v =@(lal—a,)8+(1-(3v/4nN))lal cos(C +0) , (17) 


where @(x) is the step function, @(x) = 0 for x <0 and @(x) = 1 for x >0, so 
that the tapering starts at the onset of saturation, a, = 2Y /2)73 is the 
saturation field, 6 = 4nNAy, /y, is the taper rate (6=0 for the untapered 
undulator), and y, is the resonant Lorentz factor. The modified pendulum 
equation (17) relaxes the assumption of N > 1, but still assumes y > 1 and 


allow tapering of the undulator. 
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In order to have a self-consistent theory, the dimensionless single-mode 
wave equation is also altered to account for high efficiency in terms of the 


electron phase and phase velocity [14, 15], 
a =-j<(1-(v/4nN)e“> . (18) 


In these formulations, electrons near resonance, v=0, are described 
accurately, while untrapped electrons far from resonance, v > 4nN, are 
handled less accurately. However, the untrapped electrons are expendable 
since they become randomly spread in phase as well as uncoupled from the 
bunching interaction. For electrons near resonance, the bunching rate is 
determined by the field strength lal. When lal «x, the optical field is 
considered weak, and bunching is small. When la! > 7, i.e., the optical field 
is strong, the electrons bunch rapidly, and become trapped in closed phase- 
space orbits. The current j determines the coup!ing between the electron 
beam and the optical wave. Without current, j = 0, or no bunching <..> = 0, 
there is no change in the initial optical field, a(0)=a 9. The modified 
pendulum equation (17) and wave equation (18) are valid for weak fields, 
strong fields, low gain, high gain, tapered, untapered, low efficiency, and 


high efficiency, but assume y > 1. 


In a linearly-polarized undulator, each electron experiences fast, 
periodic oscillations that are comparable to the radiation wavelength and 
modify the interaction strength. For this coupling [1], the dimensionless 
current density j should include the factor (Jo(&) - J,())? = 0.54 where 


& = K2/2(1+ K2) = 0.44 for ELF. For early ELF experiments with 450 A 
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beam current [9], the filling factor is F = 0.026, and the corrected current is 
J = 3015. 
A relativistic electron in a rectangular waveguide far above cutoff has 


its phase velocity shifted by [1] 


nNAGQ(1+K2) [p2 92 


where X is the long waveguide dimension, Y is the short waveguide 
dimension, and p and q are integers. For ELF in the TE), mode, the 
resonant phase velocity is shifted by Avo, = - 6x. A change in the undulator 
field, 5B, about the resonant field strength By also causes a change in 
phase velocity by 


4nNK? 8B 


Avp = - 
5 (1+ K?) Bo 


(20) 


Resonance in the ELF FEL is tuned using changes in the magnetic field 
5B / Bo [8}. Combining (19) and (20) shows that the waveguide increases 
the resonant undulator field strength by AB, = 240 G. 

The distribution of initial phase velocities, f(v), is important to the FEL 
growth rate and in ELF is determined by the electron beam emittance filter. 
The emittance filter results in a triangular-shaped distribution function 


with a positive slope given by 
f ,(v) = (v - v9 + 26,)/ 202 , for Vo-20,<v<vo , 


where Vo is the peak of the beam distribution. The full-width at the base of 
the trangle is cquivalent to a 6.4% energy spread so that 20, = 
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4nN x 0.064, and o, = 12 [10]. The initial input power in ELF is 30 kW in 
the TE 9, waveguide mode, and corresponds to a dimensionless microwave 


field magnitude of ay = 12. 


C. RESULTS 


Figure 5 shows the simulated gain spectrum, G(vo), and figure 6 shows 
the experimental result, G(Avg), at the end of the undulator. Recall that 
the equation (20) shows the relationship between the change in vo used in 
the simulation and the change in the magnetic field used in the experiment. 
The data is taken with J = 450 A beam current so that j = 3015, and o, = 
12 [9]. The simulation gives peak gain G,,, = 2833 at Avg = 28 away from 
resonance with a gain bandwidth of dvg = 30. The observed peak gain is 
four times higher than the calculated result, and is shifted from resonance 
more than in the simulation. However, the important qualitative features of 
the experiment are shown by the simple single-mode theory. The 
asymmetrical gain spectrum shape is in good agreement. Mode distortion 
by excitation of higher-order waveguide modes probably accounts for the 
some of higher measured gain, and the larger shift away from resonance in 
the experiment. Some of the disagreement may also be due to experimental 
error, since the powerful microwave pulse must be attenuated by many 


orders of magnitude before a measurement can be made. 
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The extraction of higher-order modes has the effect of confining the 
microwave power around the electron beam, increasing the effective filling 
factor, and increasing coupling in the later parts of the undulator. In order 
to account for mode distortion, we increase the filling factor from 0.026 to 
0.032. Figure 7 shows the simple simulation results for power evolution, 
P(t), and figure 8 shows the experiment results for power evolution, P(z). 
The later ELF experiments increase the current to 850 A [10, 11] so that 
the dimensionless current density with the new filling factor is now j = 
7000. The dimensionless microwave field is ag = 15, and o, = 12. The 
onset of saturation occurs at t= 0.5, or about 1.5 m along the undulator, 
and the efficiency is n = 5% in the untapered case. The shape of P(t) in the 
untapered case in figure 7 and figure 8 shows excellent agreement between 


the simple theory and experiment. 


In order to extend the natural saturation limit, the undulator properties 
can be altered along z to maintain resonance and interaction strength. The 
phase acceleration, 5, is given by 5 =- 4nNK2AB / B,(1+K2) when the 
undulator field strength is decreased by AB [16]. Figure 7 shows the result 
of tapering with 6 = 90x giving an efficiency of n = 31%. The experimental 
results in figure 8 shows that the efficiency is increased from 6% to 35%, 
because of tapering [11]. The efficiency is increased by about a factor of six 
in both simulation and the experiment. In both cases, tapered and 


untapered, we see excellent agreement between the experiment and the 
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simplified simulation . 


In each comparison above, the excitation of higher-order waveguide 
modes confines the microwave field tightly around the electron so that the 
filling factor F and the coupling are effectively increased. The ELF 
experiments determined that about 30% of the microwave power is in 


higher-order modes [10], consistent with this model. 


28 











690 


Figure 5. Simulated gain spectrum 
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Figure 6. Gain spectrum of the ELF experiment 
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Figure 7. Simulated gain evolution along the undulator 
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Figure 8. Gain evolution of ELF experiment 
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